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摘 要 :地 面 颤 振 模拟 试验 是 一 种 颤 振 验 证 的 全 新 试验 技术 ,可 作为 当前 赢 振 验 证 试验 手段 的 有 效 

补充 。 飞 行 器 在 实际 飞行 过 程 中 ,结构 的 动力 学 特性 及 承受 的 载荷 是 不 断 变化 的 ,对 于 受气 动 加 热 
a te 村 器 ,这 一 时 变 特 性 则 更 加 显著 。 本 研究 提出 了 基于 代理 模型 的 时 变 参 数 非 定 
气动 力 模型 建 模 方法 ,建立 了 基于 PID 控制 器 的 时 变 系 统 地 面 颤 振 试验 方法 ,并 通过 标准 试验 件 

Wa ny 试验 结果 表明 ,本 研究 提出 的 非 定常 气动 力 建 模 方法 能 够 准确 获得 具有 时 变 特 性 

的 气动 力 模型 ,建立 的 地 面 颤 振 试 验方 法 能 够 有 效应 对 颤 振 系统 的 时 变 特性 获得 准确 的 结构 颤 振 

边界 数据 。 
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Research on the ground flutter Simulation 


test method for time-varying System 
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Abstract :Cround flutter simulation test js a novel testing technology for aircraft flutter verification and an 
effective supplement to the current flutter verification test. In the process of actual flight ,aircraft’s dynamic 
characteristics and loads on the structure are constantly changing ,especially for hypersonic vehicle affected 
by aerodynamic heating. In this paper,a modeling method of time-varying parameter unsteady aerodynamic 
model is generated by the use of surrogate model. A ground flutter testing method of time-varying system 
based on PID controller is established and verified by standard component test. The testing results shows 
that the aerodynamic modeling method proposed can accurately obtain the aerodynamic model with time- 
varying characteristics, and the ground flutter testing method established can effectively respond to the 
time-varying characteristics of the flutter system and obtain accurate flutter boundary. 
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颜 振 是 弹性 结构 在 非 定常 气动 力 、 结 构 弹 性 力 
以 及 惯性 力 耦合 引发 的 气动 弹性 不 稳定 现象 ,对 于 
飞行 器 的 安全 性 具有 显著 的 影响 。 防 额 振 设计 是 飞 
行 器 结构 设计 中 的 基本 要 求 , 颜 振 边 界 获 取 也 是 飞 
行 器 强度 分 析 验 证 的 关键 环节 。 

目前 ,飞行 器 结构 颜 振 问 题 研究 中 的 试验 验证 
手段 主要 包括 风 洞 试验 与 飞行 试验 。 由 于 风 洞 试验 
通常 采用 缩 比 模型 进行 测试 , 颜 振 验证 流程 中 缺失 
了 全 机 地 面试 验 的 环节 ,造成 全 机 颤 振 飞行 试验 科 

目 多 、 周 期 长 .风险 高 。 地 面 颤 振 模 拟 试验 技术 作为 
近年 来 发 展 的 一 种 以 真实 结构 作为 验证 对 象 的 半 物 
理 仿真 验证 技术 ,具有 风险 小 、 周 期 短 、 精 度 高 的 优 
势 ,能 够 有 效 填补 目前 额 振 试验 验证 流程 中 缺失 的 
一 环 ,具有 很 大 的 工程 价值 。 

二 Keams 率 先 提 出 了 颜 振 地 面 模拟 试验 的 概 
您 ,并 对 试验 方法 与 理论 进行 了 初步 研究 。 近 年 来 ， 
国 骨 外 针对 地 面额 振 模拟 试验 中 的 关键 技术 开展 了 
诸多 研究 。Zeng 等 根据 频 域 非 定常 气动 力 计算 
方法 建立 了 试验 可 用 的 气动 力 模型 ,并 通过 重 棒 控 
制 器 实现 了 气动 力 的 精准 加 载 。Wu 等 "采用 前 置 
反馈 控制 方法 实现 了 多 激 振 器 耦合 系统 的 解 耦 控 
短 半 王 彬 文 证 提 出 了 基于 CFD 降 阶 模型 的 气动 力 
建 醒 方法 。 但 是 目前 地 面额 振 模 拟 试验 技术 的 相关 
研究 仍然 基于 系统 定常 假设 ,而 飞行 器 在 实际 飞行 
过 程 中 ,诸如 气动 加 热 对 结构 模 态 特性 的 影响 或 者 
飞 季 器 飞行 环境 变化 均 会 使 额 振 系统 表现 出 显著 的 
时 万 特 性 , 从 而 使 部 分 原 有 试验 方法 不 再 适用 中。 
基 王 上 述 原因 ,本 研究 希望 通过 开展 关于 地 面额 振 
模 秩 试验 中 的 时 变 问 题 处 理 方法 的 研究 ,包括 变 参 
数 非 定 常 气动 力 建 模 方法 以 及 颜 振 边 界 跟踪 方法 ， 
为 地 面 颤 振 模 拟 试验 提供 方法 储备 ,拓宽 地 面额 振 
模拟 试验 的 应 用 范围 。 


1 地 面 颤 振 模 拟 试验 原理 


飞行 器 地 面 闸 振 模拟 试验 是 一 种 半 物 理 仿真 试 
验 系 统 ,其 核心 思想 是 采用 激 振 设备 模拟 结构 在 飞 
行 过 程 中 承受 的 非 定常 气动 力 从 而 达到 颤 振 边界 测 
试 的 目的 。 试 验 利用 传感器 测量 结构 的 振动 响应 信 
号 ,将 其 代入 预先 计算 得 到 的 非 定 常 气动 力 降 阶 模 
型 获得 结构 的 等 效 气 动力 ,再 采用 激 振 才 将 等 效 气 
动力 实时 加 载 在 试验 结构 上 。 通 过 对 试验 结构 施加 
初始 扰动 并 观察 其 响应 , 即 可 实现 对 真实 结构 颤 振 


界 的 测试 。 试 验 基 本 原理 如 图 1 所 示 。 
计算 系统 


央 


气动 力 模型 


] | 控制 器 


功率 放大 器 


图 1 地 面额 振 模拟 试验 原理 氏 


Fig.1 Ground flutter simulation test schematic diagram 


2 时 变 系统 非 定常 气动 力 建 模 


在 地 面 颤 振 模 拟 试验 系统 中 , 非 定常 气动 力 模 
型 需要 根据 测量 得 到 的 结构 振动 响应 计算 结构 受到 
的 非 定常 气动 力 ,并 将 气动 力 减 缩 集 中 到 数 个 激 振 
点 位 置 ,测量 得 到 的 结构 振动 响应 是 其 输入 ,各 个 减 
缩 集中 的 激 振 点 气动 力 是 其 输出 。 为 了 确保 气动 力 
计算 与 加 载 的 实时 性 ,目前 研究 均 采 用 升力 面 理论 
构建 结构 与 气动 力 的 状态 空间 方程 来 实现 试验 气动 
力 建 模 ,此 时 非 定常 气动 力 的 计算 被 简化 为 数 个 低 
阶 和 矩阵 的 简单 运算 。 

虽然 在 常规 颜 振 分 析 中 气动 力 模 型 不 受 结构 动 
力学 特性 影响 ,但 地 面额 振 模拟 试验 采用 的 气动 力 
模型 需要 通过 结构 插值 矩阵 进行 降 阶 处 理 , 此 时 如 
果 飞 行 器 的 结构 模 态 特性 或 所 处 气动 环境 出 现时 变 
特性 , 则 采用 和 常 值 气动 力矩 阵 的 力 模型 将 不 再 适用 。 
因此 ,如 何 使 非 定常 气动 力 模型 在 保证 计算 效率 的 
同时 ,能 够 准确 反映 时 变 系统 与 非 定常 气动 力 的 作 
用 关系 ,是 时 变 系统 地 面 颤 振 模拟 试验 需要 重点 研 
究 的 内 容 。 


2.1 非 定常 气动 力 建 模 方法 
常规 地 面 闸 振 试验 中 的 非 定常 气动 力 建 模 方法 
已 在 文献 [6] 中 有 过 详细 推导 ,本 研究 仅 作 简要 


介绍 。 
对 于 激 振 点 / 拾 振 点 缩聚 后 的 被 测 结构 而 言 ,可 
用 降 阶 后 的 气动 力 影 响 系 数 矩 阵 4( 如 描述 试验 选 
取 的 结构 节点 位 移 zx 与 其 承受 气动 力 人 的 关系 , 具 
体 如 式 (1) 所 示 。 
f=q,A(k)z, (1) 


一 


式 中 :9。 为 来 流动 压 ;# 为 减 缩 频 率 , 旦 =wb/V。 参 
考 不 同 的 升力 面 理论 (如 偶 极 子 格 网 法 \ 活 塞 理论 )， 
即 可 构建 出 相对 应 的 气动 力 影响 系数 和 矩阵 4 (有 ) 。 
为 实现 时 域 仿真 ,需要 将 式 (1) 中 的 频 域 气动 
影响 系数 矩阵 A4(%) 转化 至 时 域 。 采 用 最 小 状态 
法 得 到 A(k) 的 有 理 函 数 近似 表达 式 ,代入 式 (1) 并 
进行 拉 氏 反 变换 ,可 得 


2 
f (2) =4, (4 二 5 + 


a (2) 
gD (7 了 | Ei. 

图 2 为 气动 力 模型 ( 式 2) 的 物理 实现 。f.() 为 
时 域 方程 ,输入 当前 时 刻 各 个 节点 的 位 移 、 速 度 与 加 
速度 即 可 求 得 当前 时 刻 的 节点 气动 力 。 在 计算 中 ， 
AA1 .4,.D\E 和 均 为 事先 求解 得 到 的 气动 力 影 
响 闲 数 矩 阵 ,经 过 激 振 点 / 拾 振 点 缩聚 ,上 述 和 矩阵 维 
度 通常 不 大 ,因此 能 够 实现 地 面 颜 振 模拟 试验 要 求 
的 时 域 气动 力 实 时 计算 和 加 载 的 要 求 。 


2 气动 力求 解 流 程 


Es Fig.2 Aerodynamic force solving process 


2 七 时 变 气动 力 影响 系数 矩阵 构建 


由 2.1 节 可 知 ,在 常规 结构 地 面 颤 振 模拟 试验 
中 气动 力 模型 的 影响 系数 和 矩 阵 是 确定 的 ,但 是 如 果 
考虑 结构 或 者 气动 力 的 时 变 特 性 , 即 结构 的 模 态 特 
性 或 所 处 气动 环境 参数 随时 间 的 改变 均 会 使 得 颤 振 
发 生 时 的 减 缩 频 率 发 生 改 变 ,进而 导致 降 阶 后 的 气 
动力 模型 发 生变 化 。 基 于 上 述 原因 ,时 变 系 统 的 地 
面 闸 振 模拟 试验 无 法 像 定常 地 面额 振 模 拟 试验 一 
样 ,采用 固定 的 矩阵 参数 进行 气动 力 建 模 。 

从 提升 试验 精度 的 角度 出 发 ,系统 在 试验 进程 
的 每 一 个 时 间 步 长 ,都 应 采用 2. 1 节 所 述 方法 求解 
一 次 当前 时 间 节 点 下 的 相关 和 矩阵 参数 ,从 而 建立 时 
变 参 数 气动 力 模型 ,但 考虑 到 计算 效率 ,这 种 方法 显 
然 无 法 满足 地 面 颤 振 模 拟 试验 系统 对 气动 力 加 载 实 
时 性 的 要 求 。 为 了 兼顾 气动 力 模型 的 精度 与 效率 ， 
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本 研究 提出 了 一 种 基于 代理 模型 的 变 参数 非 定 常 气 
动力 建 模 方法 。 

代理 模型 (surrogate model)' 是 指 对 设计 空间 
进行 抽样 后 建立 的 计算 量 远 小 于 原 模 型 而 计算 精度 
与 原 模型 相当 的 一 种 近似 数学 模型 。 代 理 模 型 建立 
后 可 以 在 分 析 过 程 中 代替 原 模型 ,从 而 解决 分 析 中 
直接 调用 高 可 信和 度 分 析 模 型 所 面临 的 计算 量 过 大 、 
计算 效率 过 低 的 问题 。 对 于 本 研究 而 言 , 原 有 高 可 
信和 度 模型 即 为 2.1 节 所 述 的 气动 力 影 响 系 数 矩 阵 计 
算 模 型 ,通过 建立 代理 模型 (此 时 气动 力 影 响 系 数 算 
阵 求解 过 程 简 化 为 一 系列 参数 已 知 的 函数 计算 ) , 即 
可 实现 在 试验 进程 的 每 一 个 时 间 步 长 中 都 求解 一 次 
当前 时 间 节 点 下 的 相关 和 抢 阵 参数 , 满足 地 面 闸 振 模 
拟 试 验 对 计算 效率 与 精度 的 需求 。 

在 航空 航天 领域 , Kriging 模型 是 最 具 代 表 性 的 
代理 模型 之 一 。 该 模型 具有 良好 的 非 线 性 近似 能 
及 多 峰值 问题 处 理 能 力 , 对 于 本 研究 需要 处 理 的 气 
动力 影响 系数 矩阵 这 类 弱 非 线性 系统 而 言 ,只 要 合 
理 选取 样本 点 即 可 获得 较 好 的 插值 近似 效果 。 

Kriging 模型 是 一 种 基于 随机 过 程 的 均 方 误差 
最 小 的 无 偏 估计 模型 。 以 任 一 时 变 气 动力 影响 系数 
和 矩阵 元 素 4;(b) 为 例 , 在 任意 时 刻 刀 其 对 应 的 值 被 随 
机 函数 4;(t) 代替 ,4;(b 只 是 该 随机 函数 可 能 的 结 
果 之 一 , 即 


LD) = PBLD) +Z00) (3) 


全 


式 中 :f,(?) 为 基 范 数 ; > Bf,(t) 为 4;(1) 的 数学 期 


望 ;Z(t) 是 均值 为 0, 方 差 为 o; 的 静态 随机 过 程 ,在 
定义 的 时 间 范 围 内 任意 2 个 取样 时 刻 的 随机 量 
Z(4) 与 Z(4) 的 协 方差 可 以 表示 为 
Cov[Z(4),2(1)] =03R(t,b) (4) 
式 中 R(t ,bs ) 为 相关 函数 。 该 函数 只 与 时 间 点 4 和 
i 之 间 的 距离 有 关 , 且 满足 : 当 距 离 趋 近 于 0 时 ,天 
数值 趋 近 于 1; 当 距 离 增加 时 ,函数 值 光滑 减 小 ; 当 
距离 趋 于 无 穷 时 ,函数 值 趋 近 于 0; 至 少 一 阶 可 导 。 
Kriging 模型 常用 的 相关 函数 包括 高 斯 函数 、 三 次 样 
条 函数 等 多 种 类 型 ,相关 函数 的 选择 可 直接 影响 算 
例 的 建 模 精度 ” 。 
对 于 本 研究 处 理 的 一 维 变量 模型 ,根据 式 (3) 
可 将 nn 个 样本 点 的 响应 定义 为 
[4;(i)], = [ 45 ) A;(i,) A;(1,)] 
=BoF +z (5) 
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式 中 
F=[f(1) fl1) … f(1,)], 
z=[Z(1) 2Z(4,) 2Z(1,)] 
根据 相关 函数 R(4 ,ts) ,可 定义 表示 n 个 样本 
点 之 间 相 关 性 的 相关 和 矩阵 R 以 及 表示 任意 一 个 样本 点 
与 其 他 所 有 样本 点 相关 性 的 相关 向 量 r, ,分 别 如 下 。 
R= (R(t,t)), eR"”" (6) 
r=[R(i,t) R(t,,t') R(t,,t')] (7) 
式 中 ,teltl 1 1 tb,]o 
样本 点 之 外 的 时 间 节 点 函数 值 可 由 已 知 样本 点 
函数 值 的 线性 加 权 得 到 , 即 
Aj(t) =cT [4,(t)], (8) 
为 满足 Kriging 模型 无 偏 估 计 E(2Z(t)) =0 的 要 
求 5 根据 式 (5 ) 与 式 (8) 可 得 
CQ rECc? 04(D],) =cfFB 
二 | E(A(D)) = [A "8 
根据 式 (3) . 式 (8) 与 式 (9), 可 得 Kriging 模型 
的 预测 值 误差 为 
Aj(t) -di(D=crz -Z(t) (10) 
C9 根据 式 (10) ,可 得 Kriging 模型 预测 值 的 均 方 误 


Fc,=f(t) (9) 


CNMSEI4(D] =E[ (4;(t) -A,(t))’] 
之 =o2(1+cTRc -2cT7r) 
因此 ,根据 Kriging 模型 预测 值 满足 无 偏 估 计 且 
9 坟 误差 最 小 的 要 求 ,该 代理 模型 的 构建 即 为 在 
FSe, =/(1) 的 前 提 下 ,寻找 满足 均 方 误差 最 小 的 c,， 
演 王 引入 拉 格 朗 日 乘 数 构建 及 (c,, 和 A) 如 下 。 
Fe ,MA)=[o(l+erRe -2cTrr)]- 
ALFTc -f(t)] 
五 分 别 对 c, 和 和 求 偏 导数 ,可 解 得 满足 要 求 的 
c, 为 
c,=R (rr,-F (FRF) (FIR, -f(t))) (13) 
显然 ,由 前 文 对 式 (13) 中 相关 参数 的 定义 可 
知 ,c, 为 阶 数 与 样本 点 数量 n 相同 的 列 向 量 ,在 此 记 
为 c = Cxi(7,)。 
将 式 (13) 代 入 式 (8) ,可 得 Kriging 模型 在 任意 
时 刻 ;的 预测 值 为 
Aj(1) =Ch (7 ) [A,()], (14) 
由 式 (5) 与 式 (13) 可 知 , 式 (14) 中 除 r, 外 其 他 
参数 只 与 样本 点 数据 有 关 , 因 此 可 以 提前 计算 得 到 ， 
而 r, 需要 将 当前 时 刻 上代 人 相关 函数 R(,,t) 求 解 ， 


(11) 


(12) 
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模型 的 输入 参数 为 时 间 , 因 此 基于 该 代理 模型 可 以 
获得 任意 时 刻 的 气动 力 。 


3 ”时 变 颤 振 边界 跟踪 测试 方法 


在 时 不 变 系统 地 面额 振 试验 中 ,通过 不 断 调 整 
风速 参数 ,观测 不 同 风速 下 的 系统 响应 , 即 可 确定 结 
构 的 闸 振 临界 速度 。 但 是 ,如 果 将 时 变 特性 引入 地 
面 颤 振 试验 系统 中 ,此 时 结构 的 颤 振 边界 也 具有 时 
变 特性 , 原 获取 颤 振 边 界 的 方式 将 不 再 适用 。 在 该 
情况 下 一 种 折 中 的 试验 方法 是 采用 某 一 速度 (或 速 
度 函 数 ) 进行 试验 ,从 而 验证 在 已 知 的 飞行 包 线 下 结 
构 不 会 发 生 颤 振 ,但 这 种 试验 方法 仅 能 作为 一 种 定 
性 的 验证 手段 ,作为 闸 振 试验 却 无 法 获得 被 测 对 象 
的 颤 振 临界 速度 ,因而 无 法 满足 工程 要 求 。 

基于 上 述 原因 ,本 研究 提出 一 种 基于 PID 控制 
的 地 面 颤 振 试验 颤 振 边界 跟踪 方法 。 该 方法 充分 利 
用 了 地 面 颤 振 试验 高 安全 性 的 优势 ,依据 结构 响应 
变化 趋势 调整 风速 ,能 够 较为 准确 地 得 到 结构 的 闸 
振 临 界 速 度 。 

在 地 面 闸 振 试验 中 ,人 为 调整 风速 参数 寻找 被 
测 结构 颤 振 边界 的 逻辑 可 以 概括 为 :结构 响应 发 散 
则 降低 风速 ,结构 啊 应 收敛 则 增高 风速 。 如 果 将 结 
构 响 应 幅 值 变化 趋势 定义 为 A(A > 0 为 响应 发 散 ， 
入 <0 为 响应 收敛 ) , 则 显然 速度 修正 量 AV 应 与 入 保 
持 负 相 关 关 系 ,本 章 的 颜 振 边界 跟踪 系统 即 希 望 通 
过 PID 控制 器 实现 该 速度 修正 过 程 。 


3.1 响应 信号 处 理 


考虑 一 段 典 型 的 地 面 颤 振 仿真 系统 结构 响应 信 
号 A(t) 如 图 3 所 示 , 显 然 在 仿真 时 段 内 该 信号 先 发 
散 后 收敛 。 


响应 位 移 /mm 


0 5 10 Lb] 20 
时 间 /s 
图 3 结构 响应 信号 


Fig.3 Structural response signal 


为 获取 响应 幅 值 变化 趋势 ,首先 应 获得 振动 幅 
值 信号 4'(t) ,选取 的 时 间 点 # 访 
后 的 响应 幅 值 信号 如 图 4 所 示 。 
A(t,) =0 
[ia (15) 
A(t._,) >0 


应 满足 式 (15 ) ,处 理 


响应 幅 值 /mm 


| 0 3 10 13 20 
A | 时 间 / S 

CN 

[ey 到 4 阶 跃 幅 值 信号 

0) Fig.4 Step amplitude signal 


己 由 图 4 曲线 的 局 部 , 即 黑 杠 中 的 内 容 可 以 看 出 ， 
听 应 信号 起 始 阶段 虽然 幅 值 整体 变化 趋势 为 增高 ， 
全 局 部 区 域 仍 会 出 现 降低 的 变化 趋势 , 即 振动 幅 值 
信号 并 非 单调 递增 ,而 在 实际 试验 中 由 于 噪声 等 护 
动 内 素 该 现象 会 更 加 显著 。 为 了 增加 控制 系统 的 稳 
乱 隧 ,需要 排除 诸如 此 类 的 信号 局 部 扰动 影响 ,使 控 
制 器 在 最 大 程度 上 关注 幅 值 整体 的 变化 趋势 ,因此 
毕 处 对 信号 做 数据 平滑 处 理 ,采用 滑动 平均 法 加 如 
式 困 6) 所 示 ,积分 宽度 A7 为 选取 观测 值 的 宽度 ,可 
由 盟 体 算 例 决定 ,适当 增 大 该 值 可 以 增加 信号 平滑 
en 号 处 理 结果 如 图 5 所 
示 了 此 处 只 做 分 段 积分 虽然 会 导致 信号 数值 发 生变 

化 ,但 由 于 系统 仅 关注 信号 的 变化 趋势 ,因此 数值 改 
变 并 不 会 影响 控制 效果 。 
40 = 4(D)d 上 6 


# 构 响应 幅 值 整体 变化 趋势 
的 连续 信号 4(1) 。 一 个 较 小 的 时 间 间 隔 At, 则 
描述 结构 响应 ee 
(17) 表 示 。 此 处 不 直接 对 信号 4(1) 作 微分 处 理 是 
考虑 到 实际 试验 中 由 于 噪声 等 干扰 因素 会 降低 信和 号 
质量 ,直接 微分 处 理会 导致 系统 稳定 性 降低 。 

4( 

A A(i-Ai) 

采用 式 (17) 对 图 5 所 示 信 号 进行 处 理 ,所 得 结 
果 如 图 6 所 示 。 由 图 3 可 知 结构 响应 信号 先 发 散 且 


(17) 
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发 散 速度 不 断 变 慢 ,在 9 ~10s 的 区 段 内 由 发 散 变 为 
收敛 ,而 图 6 中 入 先 为 正 且 其 值 不 断 减 小 ,9 ~ 10 s 
的 区 段 内 由 正 变 为 负 并 在 之 后 一 直 为 负 , 膏 明 入 外 
够 较 好 描述 结构 响应 号 的 变化 趋势, 可 以 作为 控 
制 侨 的 输入 量 。 

1.0 


幅 值 分 段 积分 


0 5 10 15 20 
时 间 /s 


图 5 平滑 幅 值 信号 


Fig.5 Smooth amplitude signal 


.00 5 10 15 20 
时 间 /s 
了 6” 幅 值 变化 趋势 信号 


Fig.6 Amplitude variation trend signal 


3.2 控制 器 设计 
PID 控制 器 动态 方程 如 式 (18) 所 示 '”。 
sD) = Ee) + Fe dt kor! sD (18) 


对 于 PID 控制 器 设计 而 言 , 最 关键 的 步骤 就 是 

PID 参数 整定 , 即 参 考 系统 被 控 量 的 时 域 响应 调节 

式 (18) 中 的 K,、K;=K,/7T;、K,=K,7。 对 于 本 次 研究 

对 象 地 面 颤 振 试验 系统 而 言 ,调节 量 为 空 速 变 化 量 

AV(t) ,被 控 量 为 描述 振幅 变化 趋势 的 量化 参数 
A(i) ,这 二 者 间 存 在 非 线性 关系 ,可 描述 为 

A(1) =FCTCDD))+ATFCD)+aCATCD) +6(1) 

(19) 

式 中 :ogo(AV(z) ) 为 非 线性 项 ;5(b 为 干扰 项 。 由 于 

与 k( V(t) +AV(t) ) 相 比 这 两 项 均 为 小 量 ,因此 A (7) 

与 AV(t) 整体 上 保持 正 相 关 , 但 其 存在 会 导致 A(t) 
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如 图 6 所 示 的 震荡 现象 难以 消除 ,从 物理 意义 上 讲 
只 要 A(t) 保 持 在 0 附近 较 小 的 区 间 内 都 可 以 视 系 
统 响 应 信号 为 等 幅 振 动 ,此 时 和 (i) 本 身 的 震荡 并 没 
有 实际 意义 。 因 此 ,在 PID 参数 整定 的 过 程 中 应 选 
取 AV(zt) 作 用 后 的 空 速 V(z) 而 非 和 A(z) 作为 参考 信 
号 ,该 信号 同时 也 是 地 面额 振 试验 中 希望 得 到 的 闸 
振 边 界 预测 结果 信号 。 


4 ”地 面 颤 振 试验 仿真 平台 构建 


为 验证 前 文 所 述 变 参数 非 定常 气动 力 模型 及 字 
振 边界 预测 方法 的 可 行 性 ,本 章 在 matlab 中 的 simu- 
link 模块 下 搭建 试验 系统 的 仿真 平台 , 并 采用 一 模 
型 算 例 对 仿真 平台 进行 验证 。 


广 模型 介绍 
CO 本 仿真 算 例 及 后 续 试验 均 选用 一 典型 亚 声速 机 


af 
惨 际 准 试验 件 为 研究 对 象 , 该 模型 由 金属 梁 架 及 维 形 


术 框 组 成 ,其 中 承 力 梁 架 材料 为 7075 铝 合金 ,约束 方 


纹 泗 根 固 支 , 模 型 的 主要 模 态 参数 如 表 1 所 示 。 
表 1 机 杜 前 三 阶 模 态 频率 


Lower three natural frequencies of the wing 


Tab. 1 


异 态 阶 数 模 态 名 称 模 态 频率 /Hz 
1 一 阶 弯 曲 3.44 
2 一 阶 扭转 9.09 
3 二 阶 弯曲 11.92 


过 赢 振 分 析 


所 频 域 颜 振 分 析 采 用 MSC. Nastran 求解 。 为 了 模 


拟 气 动力 的 时 变 特性 ,研究 通过 修改 来 流 密度 实现 ， 
选取 来 流 密度 比 由 0.9 线性 降 至 0.7 ,在 区 间 内 均匀 
选取 5 个 采样 点 进行 闸 振 分 析 。 图 7 为 密度 比 为 
0.7 时 的 vg 与 v7 图 ,得 到 各 个 采样 点 颤 振 临 界 速 
度 及 颤 振 频 率 如 表 2 所 示 。 
表 2 采样 工 况 颤 振 特性 ( 有 限 元 ) 
Tab.2 Sampling point flutter characteristics ( FEM) 


密度 比 颤 振 速度 /(m .ss') 颤 振 频率 /Hz 
0.90 50.0 | 
0. 85 S12 7.10 
0. 80 52.6 7.10 
0.75 54.1 7.09 
0.70 55.7 7.09 


Pei As 
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-一 阶 模 态 
-二 阶 模 态 
-下 三 阶 模 态 


20 30 40 50 60 70 80 
速度 wm 。s 7) 
(a) v-g 曲线 


频率 1/Hz 


2 
20 30 40 50 60 70 80 


CD) wy 曲线 


图 7 vw-g 与 vf 曲线 图 (密度 比 =0.7) 
Fig.7 wvw-g and v-f plots( density ratio =0.7) 


4.3 平台 构建 


采用 本 研究 第 2、3 节 方 法 ,参考 图 2 构建 气动 
力 仿真 平台 如 图 8 所 示 。 其 中 ,具有 时 变 特性 的 气 
动力 影响 系数 矩阵 由 子 模块 计算 得 到 , 子 模块 中 包 
括 了 自 变 量 为 时 间 的 Kriging 模型 运算 ( 式 14) 及 考 
虑 来 流速 度 的 矩阵 系数 运算 ; 速度 控制 器 通过 读 取 
输入 气动 力 模型 的 结构 响应 计算 修正 风速 ,并 将 计 
算 结果 输出 给 各 个 气动 力 影 响 系 数 和 矩阵 子 模块 , 控 
制 器 可 由 开关 控制 是 否 将 其 接 入 仿真 系统 ,在 不 接 
入 则 以 固定 风速 值 进行 仿真 。 

与 实际 的 地 面 颤 振 模拟 试验 相 比 ,仿真 系统 没 
有 试验 对 象 提供 结构 响应 信号 ,因此 还 需要 建立 结 
构 响 应 模块 使 仿真 系统 形成 闭环 。 通 过 模 态 参数 建 
立 结 构 的 状态 空间 方程 ,结构 响应 模块 构建 如 图 9 
所 示 , 其 中 通过 “phi_f 与 "phi_z 两 步 矩 阵 运算 实现 
气动 力 模 型 与 结构 状态 空间 方程 间 物 理 坐 标 与 广义 
坐标 的 转换 。 


Outl 
dynamic pressure 


7 x'=Ax+Bi 
AL y=Cx+Du 

1 人 
Gy phif structrual dynamics 


CO 图 9 结构 响应 计算 模块 


OO Fig.9 Stuctural response calculation module 


phi z 


42 仿真 验证 


和 ?构建 仿真 平台 的 气动 力 模型 时 设 定 每 两 个 采样 
站 
籽 0~200 s。 

_ 伍 首先 将 仿真 平台 子 模块 中 的 时 间 输 出 设 为 固定 
信和 郝 进行 定常 地 面 颤 振 仿真 测量 各 采样 工 况 下 的 
赢 所 临界 速度 ,通过 不 断 调 整风 速 观 测 结构 模块 输 
出 网 响 应 ,得 到 仿真 结果 并 与 有 限 元 数值 计算 结果 
( 案 ? ) 对 比 误差 如 表 3 所 示 , 可 知 该 颤 振 仿 真 系统 
在 0 ~200s 的 过 程 中 颤 振 临界 速度 由 51.4 m/s 升 
至 56.9 m/s。 在 理想 情况 下 ,采用 该 仿真 平台 所 得 
的 颤 振 临界 速度 应 与 有 限 元 数值 计算 结果 保持 一 
致 ,但 4.3 节 所 述 构建 结构 响应 计算 模块 的 状态 空 
间 建 模 方法 为 降低 计算 规模 会 忽略 高 阶 模 态 ,从 而 
导致 了 表 3 所 示 的 仿真 误差 。 

表 3 采样 工 况 颤 振 特性 (仿真 ) 


Tab.3 Sampling point flutter characteristics (simulation) 


二 仿真 结果 / 数值 计算 结果 / 误差 / 
密度 比 | | 
(ms  ) (ms  ) % 
0.90 51.4 50.0 2.8 
0. 85 52.3 51.2 2.5 
0.80 53.8 52.6 2.3 
0.75 33.3 54.1 2.2 
0.70 56.9 55.7 2.2 
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人 vel_control 


fix_vel 


图 8 气动 力 计算 模块 


Fig.8 Aerodynamic force calculation module 


再 将 仿真 时 间作 为 输入 信号 引入 系统 子 模块 ， 
速度 控制 器 在 仿真 开始 10 s 后 接 入 系统 ,其 初始 输 
出 风速 设 定 为 ”= 52 m/s。 仿 真 所 得 结构 响应 以 及 
跟踪 速度 曲线 如 图 10 所 示 ,结合 表 3 定常 仿真 数据 
(图 10 中 三 角形 标记 点 ) 可 以 看 出 :在 仿真 前 10s 由 
于 速度 控制 器 尚未 接 和 人 而 初始 风速 大 于 系统 此 时 的 
颤 振 临 界 速度 ,结构 响应 快速 发 散 ; 速度 控 制 器 在 
10 s 接 入 系统 后 对 信和 号 的 发 散 趋 势 作出 反应 ,输出 
风速 开始 降低 ,使 得 结构 响应 在 10 ~20 s 间 由 发 散 
转 为 收敛 ;在 30 ~ 200 s 的 区 间 内 ,结构 响应 信号 只 
有 幅度 较 小 的 缓慢 收敛 或 发 散 , 整 体 上 响应 幅 值 变 
化 不 大 ,从 额 振 分 析 的 角度 出 发 可 以 将 此 时 段 的 结 
构 视 为 闸 振 临界 状态 , 而 这 一 区 段 内 速度 控制 器 输 
出 的 速度 也 与 数值 计算 的 颤 振 临 界 速 度 基本 一 致 ,说 
明 该 仿真 平台 能 够 表现 出 颤 振 系 统 的 时 变 特性 且 速 
度 追踪 系统 能 够 得 到 系统 颜 振 临界 速度 的 时 变 历程 。 
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图 10 ”时 变 地 面 颤 振 系统 仿真 结 


Fig. 10 Time-varying ground flutter system simulation results 
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综 上 所 述 ,在 仿真 分 析 中 气动 力 模型 能 够 较 好 
表现 其 时 变 特性 ,速度 控制 融 也 能 够 实现 第 3 节 所 
述 的 参数 调节 过 程 ,跟踪 效果 良好 , 且 仿真 平台 整体 
运行 速度 较 快 ,能 够 满足 地 面 颤 振 模 拟 试验 系统 对 
于 计算 实时 性 的 要 求 。 


5 地面 颤 振 验证 试验 


5.1 试验 系统 集成 


为 了 便于 地 面额 振 试 验 的 物理 实现 ,试验 前 首 
先 需要 选取 合理 的 激 振 点 / 拾 振 点 位 置 ,以 便利 用 少 
量 的 激 振 点 和 拾 振 点 反应 整体 结构 气动 特性 。 本 研 
究 设 定 激 振 点 / 拾 振 点 数目 为 4, 采 用 基于 广义 力 等 
效 的 遗传 算法 进行 优化 配置 中 ,为 了 将 激 振 点 布置 

站 承 力 结构 上 ,将 优化 的 位 置 在 框架 上 作为 约束 
依从 ,同时 为 了 满足 时 变 系统 整个 时 段 的 模拟 精度 
要 惑 ,将 各 个 时 段 的 广义 力 等 效 目 标的 加 权 作 为 束 
体 优化 目标 , 进行 综合 优化 ,最 终 优 化 激 振 点 位 置 
(名 点 与 激 振 点 共 点 ) 如 图 11 所 示 。 
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图 11 激 振 点 位 置 
Fig.11 Location of excitation point 

在 每 一 个 优化 所 得 节点 处 安装 激 振 设备 及 数据 
采集 设备 ,构建 试验 系统 如 图 12 所 示 。 试 验 采 用 激 
光 位 移 传感器 测 得 节点 位 移 并 将 信号 微分 得 到 节点 
速度 ,采用 加 速度 传感器 测 得 节点 加 速度 ;在 激 振 器 
顶 杆 末端 安装 力 传感器 测量 实际 加 载 的 力 信和 号, 通 
过 重 棒 控制 进行 激 振 力 的 实时 反馈 控制 ,确保 激 

振 器 能 够 精确 模拟 等 效 气动 力 。 


5.2 试验 验证 


验证 试验 首先 采用 固定 风速 参数 进行 ,验证 时 
变 非 定常 气动 力 模型 实际 效果 。 试 验 参 数 与 4.4 市 
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仿真 系统 一 致 ,每 两 个 采样 点 间 的 时 间 间 隔 为 50 s， 
试验 及 采样 周期 均 为 200 s, 获 得 的 监测 点 位 移 如 图 
13 所 示 。 
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图 12 ”地面 颤 振 模拟 试验 系统 
Fig. 12 Ground flutter simulation test system 
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图 13 固定 速度 下 的 试验 结果 
Fig. 13 Test results with fixed velocity 

由 试验 结果 可 知 ,在 v=52 m/s 的 固定 风速 下 ， 
结构 响应 首先 发 散 ,在 20 s 前 后 转 为 收敛 且 收敛 速 
度 逐 渐变 快 ,这 与 理论 计算 结果 推导 出 的 响应 变化 
趋势 基本 一 致 ,表明 基于 代理 模型 的 时 变 参数 非 定 
常 气动 力 建 模 方法 能 够 充分 体现 系统 的 时 变 特征 ， 
可 用 于 地 面 颤 振 模 拟 试验 。 

为 验证 颤 振 边界 跟踪 算法 ,采用 相同 的 非 定常 
气动 力 模型 ,引入 颤 振 边 界 跟踪 模块 重新 进行 试验 ， 
验证 颤 振 边界 跟踪 方法 的 效果 。 在 相同 的 试验 设置 
条 件 下 进行 试验 ,同样 以 v =52 m/s 的 初始 空 速 开展 
试验 ,得 到 结构 监测 点 位 移 与 颤 振 临 界 速 度 曲 线 如 图 
14 所 示 ,其 中 颤 振 临 界 速度 曲线 图 中 实 线 为 试验 所 得 
数据 ,虚线 为 图 10 中 仿真 系统 速度 跟踪 所 得 数据 。 


响应 位 移 /mm 


0 50 100 150 200 
时 间 /s 
60 
一 一 试验 数据 
ss 仿真 数据 一 一 一 一 一 


速度 /m。s 


0 50 100 150 200 
时 间 /s 
图 14 采用 速度 跟踪 系统 的 试验 结果 


Fig.14 Test results with the usage of velocity tracking system 

如 图 14 所 示 , 试 验 初期 结构 响应 快速 发 散 , 在 
试验 开始 10 s 时 开启 赢 振 边界 跟踪 模块 ,可 观测 到 
结构 响应 由 缓慢 发 散 逐 步 变 为 等 幅 振 沪 , 阁 振 临界 
i 试验 结构 


嘛 如 与 仿真 结果 (图 10) 相 比 ,并 未 出 现在 开启 速度 
探 重 器 后 结构 响应 由 发 散 转 变 为 收敛 的 过 程 , 取 而 
代 沁 的 是 较为 平缓 的 幅 值 变 化 趋势 ,经 重复 试验 与 
企 辽 可 知 该 区 别 的 原因 在 于 被 测 结构 由 于 组 成 材料 
及 向 型 复杂 等 因素 存在 较 强 非 线性 ,因此 与 仿真 所 
2 


缚 丽 对 于 风速 的 变化 并 不 敏感 。 
二 此 外 由 图 14 中 速度 跟踪 曲线 可 知 , 该 颜 振 边界 
追踪 方法 在 理想 模型 仿真 与 实际 试验 中 均 能 够 较 好 
的 咏 踪 颤 振 系 统 的 临界 速度 , 可 用 于 具有 显著 时 变 
特 轰 的 地 面 颜 振 模拟 试验 系统 。 

-人 选取 采样 点 所 处 时 间 的 速度 追踪 数据 (50 s、 
10 反 .150 s.200 s) 并 与 有 限 元 数值 计算 结果 对 比如 
表 4 所 示 。 由 于 颤 振 数值 分 析 方法 难以 准确 考虑 真 
实 结构 存在 的 阻尼 特性 ,试验 结果 整体 略 大 于 数值 
计算 结果 是 符合 预期 的 ,该 现象 也 印证 了 采用 真实 
结构 进行 的 地 面 颤 振 模 拟 试验 在 颤 振 研究 领域 的 独 
特 优势 。 


表 4 采样 工 况 颤 振 特 性 (试验 ) 


Tab.4 Sampling point flutter characteristics (test) 


采样 时 间 点 / 试验 结果 / ”数值 计算 结果 / ”误差 / 
密度 比 
S (ms ) (ms ) % 
50 0.85 52:9 $1:2 3:2 
100 0.80 54.4 52.6 3.4 
150 0.75 56.0 54.1 3.4 
200 0.70 S75 55;7 3;3 
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6 结 论 


1) 提 出 基于 代理 模型 的 时 变 参数 非 定常 气动 力 
建 模 方 法 ,并 通过 试验 验证 可 充分 描述 颤 振 系 统 中 
气动 力 模型 的 时 变 特性 。 

2) 建立 了 采用 PID 控制 器 的 颤 振 边界 预测 方 
法 ,并 推导 了 相对 应 的 信号 处 理 方法 。 试 验证 明 ,该 
方法 能 够 实时 读 取 时 变 地 面 闸 振 试验 中 研究 对 象 在 
当前 时 刻 的 颜 振 临 界 速度 ,获取 整个 时 变 历程 中 试 
验 对 象 的 颤 振 临 界 速度 变化 曲线 ,追踪 精度 较 高 ,可 
用 于 具有 显著 时 变 特性 的 地 面 颤 振 模 拟 试验 系 统 ， 
弥补 当前 试验 方法 的 缺陷 。 

3) 时 变 系统 地 面 闸 振 试验 测试 结果 与 数值 计算 
结果 一 致 性 较 好 ,试验 系统 运行 效率 较 高 ,证 明了 本 
研究 提出 的 时 变 气 动力 模型 和 赢 振 边界 跟踪 算法 适 
用 于 时 变 系统 的 闸 振 测试 。 
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